The geometry of a naked singularity created by standing waves near a Schwarzschild 
horizon, and its application to the binary black hole problem 
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The most promising way to compute the gravitational waves emitted by binary black holes (BBHs) 
in their last dozen orbits, where post-Newtonian techniques fail, is a quasistationary approximation 
introduced by Detweiler and being pursued by Price and others. In this approximation the outgoing 
gravitational waves at infinity and downgoing gravitational waves at the holes' horizons are replaced 
by standing waves so as to guarantee that the spacetime has a helical Killing vector field. Because 
the horizon generators will not, in general, be tidally locked to the holes' orbital motion, the standing 
waves will destroy the horizons, converting the black holes into naked singularities that resemble 
black holes down to near the horizon radius. This paper uses a spherically symmetric, scalar-field 
model problem to explore in detail the following BBH issues: (i) The destruction of a horizon 
by the standing waves, (ii) The accuracy with which the resulting naked singularity resembles a 
black hole, (iii) The conversion of the standing-wave spacetime (with a destroyed horizon) into a 
spacetime with downgoing waves by the addition of a "radiation-reaction field" . (iv) The accuracy 
with which the resulting downgoing waves agree with the downgoing waves of a true black-hole 
spacetime (with horizon) . The model problem used to study these issues consists of a Schwarzschild 
black hole endowed with spherical standing waves of a scalar field, whose wave frequency and near- 
horizon energy density are chosen to match those of the standing gravitational waves of the BBH 
quasistationary approximation. It is found that the spacetime metric of the singular, standing-wave 
spacetime, and its radiation-reaction-field-constructed downgoing waves are quite close to those for a 
Schwarzschild black hole with downgoing waves — sufficiently close to make the BBH quasistationary 
approximation look promising for non-tidally-locked black holes. 

PACS numbers: 04.25.-g, 04.25.Nx, 04.30.-w 



I. INTRODUCTION AND SUMMARY 

It is very important, in gravitational astronomy, to 
have accurate computations of the gravitational waves 
from the inspiral of a black hole binary [ij. However, 
computing these waves is extremely challenging: for the 
last w 25 cycles of inspiral waves, post-Newtonian ap- 
proximations fail 0], and numerical relativity can not 
yet evolve the full dynamical equations in this regime. It 
appears that the best hope for accurately computing the 
BBH inspiral waves is by a quasi-stationary approxima- 
tion In this approximation, the energy and angular 
momentum of the binary are conserved by the imposition 
of standing gravitational waves, and the spacetime has a 
helical Killing vector field. The standing-wave radiation 
required to keep the orbit stationary is computed by de- 
manding that the energy contents of the gravitational 
waves (GW) be minimized Q. 

Standing-wave radiation consists of a sum of ingoing 
and outgoing radiation at infinity, and downgoing and 
upgoing radiation at the black- hole horizons. The physi- 
cal spacetime, with purely outgoing waves at infinity and 
downgoing waves at the horizons, can be recovered from 
the standing-wave spacetime by adding a perturbative 
radiation-reaction field 0. The solution for the BBH 
inspiral consists of a series of quasi-stationary solutions 
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that evolve from one to another via energy and angular 
momentum loss triggered by the radiation-reaction field. 
The waves measured at a detector can be deduced from 
this sequence of quasistationary solutions. 

The black holes comprising the binary are tidally 
locked if their horizon generators are static in the frame 
co-rotating with the orbit. In the tidally locked case, the 
metric perturbations necessary to keep the black holes 
on a stationary orbit are static in the co-rotating frame, 
and the black holes can be regarded as having bifurcate 
Killing horizons (both a past horizon and a future hori- 
zon). 

In reality, the black holes are not tidally locked. Their 
mutual tidal forces are not strong enough to maintain 
locking during the inspiral. In the absence of tidal lock- 
ing, the standing waves of the standing-wave approxi- 
mation destroy the black- hole horizons: the downgoing 
waves destroy the past horizon by building up an infi- 
nite energy density at the past horizon, and the upgoing 
waves destroy the future horizon. Therefore, we expect 
that forcing the orbit to be stationary via the addition 
of standing gravitational waves will strip the Kerr black 
holes of their horizons and leave naked singularities in 
their place |f|. 

Despite this radical change in the character of the or- 
biting bodies, it is reasonable to expect that the standing- 
wave solution will give a quite accurate approximation to 
the true physical black-hole spacetimes everywhere ex- 
cept very near the black-hole horizons. In order to verify 
or refute this expectation, it is necessary to explore the 
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nature of the singularities created by the standing gravi- 
tational waves and to test how well the physical solution 
with true black holes can be extracted from the standing- 
wave solution with naked singularities. 

As a first step in such an exploration, we consider in 
this paper a simple model problem designed to give in- 
sight into the nature of the singularities generated by 
the standing gravitational waves, and the accuracy with 
which the physical, BBH spacetime can be recovered from 
the standing- wave, singularity-endowed spacetime. 

Our model problem is a single, spherically symmetric 
black hole that is converted into a naked singularity by 
spherical standing waves of a scalar field. 

We begin our analysis in Sec.[H]by describing the map- 
ping between the BBH problem, into which we seek in- 
sight, and our spherical, scalar- field model problem. In 
particular, we deduce what should be the range of scalar- 
field amplitudes and frequencies in order to mock up the 
gravitational waves of the BBH problem. 

Then in Sec. 11111 we construct and explore the 
standing-wave spacetime for our spherical model prob- 
lem. We initially treat the standing-wave scalar field as 
residing in the unperturbed Schwarzschild spacetime of 
the black hole, and we use Regge- Wheeler first-order per- 
turbation theory to compute the scalar-energy-induced 
deviations of the hole's metric from Schwarzschild. The 
metric perturbations consist of a static component and a 
component varying in time at twice the scalar-field fre- 
quency (see Fig. |2 below). The oscillatory component is 
smaller than the static one and higher-order harmonics 
of both the field and the metric are strongly suppressed. 

The static metric perturbation grows divergently as 
one approaches the Schwarzschild horizon — an obvious 
indication of the horizon's destruction by the standing- 
wave stress-energy. To explore the structure of the re- 
sulting naked singularity, in Sec. IIII 51 we abandon per- 
turbation theory and switch to the fully nonlinear, cou- 
pled Einstein equations and scalar-field equations. To 
simplify the analysis, we focus solely on the static part 
of the singularity's metric; we do this by time averaging 
the scalar stress-energy tensor before inserting it into the 
fully nonlinear Einstein equations. We solve the resulting 
equations numerically to obtain the spacetime geometry 
outside the singularity. The geometry's embedding dia- 
gram (Fig. |3 below) and the redshift seen by a distant 
observer (Fig. 0] below) show that the spacetime remains 
nearly Schwarzschild outside the Schwarzschild horizon, 
but deviates strongly from Schwarzschild at r w 2M and 
below. (Here M is the mass of the hole-like singularity 
and we use geometrized units c = G = 1 everywhere in 
this paper.) Above r = 2M, the standing-wave space- 
time is very nearly identical to the Schwarzschild space- 
time down to radii that are well inside the inner edge of 
the effective potential (Fig.[SJ). Below r — 2M, radial dis- 
tance changes far more slowly than areal radius; i.e., g rr 
tends to as r — > 0. The redshift seen by an external ob- 
server rises rapidly when the emitter falls inside r = 2M. 
However, a signal from the singularity at r = may be 



infinitely redshifted or infinitely blueshifted, depending 
on the choice of scalar field parameters. 

In Sec. IIVI we turn to the model spherical spacetime 
that mocks up our desired BBH solution: the space- 
time of a Schwarzschild black hole with downgoing scalar 
waves. Not surprisingly, the metric perturbations in- 
duced by the downgoing scalar-wave energy are those of 
the Vaidya solution of Einstein's equations — a slowly 
growing black hole with a smooth, non-singular future 
horizon. This spacetime is well approximated, for short 
time intervals, by the Schwarzschild solution with (con- 
stant) Schwarzschild mass equal to the instantaneous 
Vaidya mass. 

Finally, in Sec.[V]we demonstrate that by adding a per- 
turbative radiation-reaction field to the standing-wave 
solution, a downgoing solution to the scalar-wave equa- 
tion can be recovered. We explore the level of agree- 
ment between these downgoing waves that live in the 
singularity-endowed standing-wave spacetime and the 
downgoing waves in the Schwarzschild approximation 
to the Vaidya spacetime. The agreement (for details 
see Sec. and Fig. below) is rather good for scalar- 
wave amplitudes and frequencies that mock up the BBH 
problem — sufficiently good to give optimism that the 
standing-wave approximation will give accurate gravita- 
tional waveforms for the final stages of binary-black-holc 
inspiral. 

II. THE MAPPING BETWEEN THE BBH 
PROBLEM AND OUR MODEL SCALAR-FIELD 
PROBLEM 

In our exploration of the quasistationary, standing- 
wave approximation for black- hole binaries we shall study 
several spherically symmetric spacetimes, each endowed 
with a standing-wave scalar field. In Sec. IIII Al the 
spacetime will be Schwarzschild, or Schwarzschild with 
first-order gravitational perturbations generated by the 
scalar- field stress-energy tensor. In Sec. IIII Bl the space- 
time will be that of a naked singularity generated by the 
coupled, time-averaged Einstein-scalar-field equations. 
In this section we shall identify the parameter regime rel- 
evant to gaining insight from these spacetimes into the 
binary black hole problem. 

In each of these spherical spacetimes, the scalar field 
must be a solution to the wave equation: 

□ $=-L(V^3 a/3< f.a),/3-0, (1) 

where g a j3 is the spacetime metric with the interval 

ds 2 = f(r,t)dt 2 + g{r,t)dr 2 + r 2 (d6 2 + sin 2 9d(f> 2 ). (2) 

We assume that the scalar field is monochromatic with 
frequency u), and we write it in the form 

* = »(*^), (3) 
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where 5ft () denotes the real part and the phase was set by 
the choice of the zero of time t. 

The scalar field $ serves as the source of curvature in 
the Einstein equations, 



G a p = 8nT, 



a/3 ■ 



(4) 



where the stress-energy tensor depends on the scalar field 
according to 

Tap = 4^,a$ /3 ~ ^9*P*,K*' M (5) 

(cf. Eq. (20.66) of 7] or Eq. (A. 11) of 0). 

We can re- write equations Q and JSJ in a simpler form 
via the Ricci tensor: 



(In the spirit of this approximate analysis we here ignore 
the gravitational blueshift of the energy.) By equating 
this energy density to the value of Too at the horizon, 
computed by inserting Eq. into Eq. 1(5} , we obtain 
the scalar-field amplitude inside the peak of the effective 
potential: 



(a/M) 



M . 



(11) 



Using equations JSJ and 1|11[1 . we can compute the de- 
sired scalar-field frequency and amplitude for the bound- 
aries of the region of interest: 



a = 6M 



0.19/M, * in w 0.017M; 



(12a) 



(6) 



Relevant ranges for the scalar-field frequency and am- 
plitude are determined by the binary black hole problem 
we are modeling. Suppose that the black holes in the bi- 
nary have equal mass M, and let a be their radial sepa- 
ration. Since we are interested in the late inspiral, where 
the post-Newtonian methods fail, the desired range of 
parameters should correspond to 6 < a/M < 15 0- 

The Keplerian orbital frequency of the black holes is 



(a/ 'Mf 



(J) 



The gravitational wave frequency is twice the Keplerian 
frequency, and we set our scalar-field frequency equal to 
the GW frequency: 



2 2 

2ft = 

MV (a Mf 



(8) 



The power going down a black hole due to the orbital 
motion of its companion is approximately 



p gw = ^mVo 6 
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(9) 



where \x is the mass of the companion 

HG3. Although 
the calculations in Refs. H E3 underlying Eq. © 

were 

carried out under the assumption jj, -C M, we will use 
Eq. I© to approximate the power for equal mass black 
holes, fi = M. This approximation is not too worri- 
some because we are interested in the general features 
of the scalar-field model, which roughly corresponds to 
the interesting range of BBH separations, rather than in 
the quantitative results for this model. We select the 
scalar-field amplitude by demanding that its energy den- 
sity near the horizon equal the GW energy density there: 



GW 



dE ^ 

dV ~ 4tt(2M) 2 



(10) 



15M 



UJ 



0.049/M, % n w 0.0011M. (12b) 



III. STANDING- WAVE SCALAR FIELD 

We now turn to the standing-wave scalar-field space- 
time that mocks up the spacetimes of the BBH standing- 
wave approximation. The metric of this spacetime has 
the form of Eq. J2J and the standing-wave scalar field 
follows from Eq. ©: 



\I/(r) cos uit 



(13) 



where ty(r) is now real. 

We shall treat the standing-wave scalar field twice 
via two different simplifying assumptions. First, in 
Sec. IIII Al we will consider the scalar field perturbatively; 
its wave equation will be that of the Schwarzschild space- 
time, and its stress-energy will generate first-order per- 
turbations of the metric away from Schwarzschild. Then 
in Sec. IIII 51 we will consider the fully nonlinear Einstein- 
scalar-field spacetime but with the scalar stress-energy 
averaged over time to make the metric static. 



A. Perturbative standing-wave solution 

1. Perturbative formalism for the standing-wave spacetime 

In our first approach, the lowest-order solution for the 
scalar field is computed by solving the wave equation Q 
in the Schwarzschild background with the metric 



ds 2 = g^pdx a dx 13 

= -(l-2/r)dt 2 + 



(14) 



1 - 2/r 



dr 2 



r 2 dQ 2 



where we rescale so that M = 1. The wave equation 
simplifies as follows (cf. Eq. (32.27b) of 0): 



d 2 ^ 

dr* 2 



-u 2 + (l-2/r)- 



(15) 
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where r* is the Regge- Wheeler tortoise coordinate ^lj, Regge- Wheeler form: 



2 In (r/2- 1) 



(16) 



Because u> 2 dominates the right hand side of Eq. i|15|) 
both far from the horizon (r 3> 2) and very near the 
horizon, the scalar field will oscillate with a nearly con- 
stant frequency lu in those regions. In between, where 
the effective potential 



V(r*) = (l-2/r)(2/r 3 ), 



(17) 



is significant, there is an intermediate transitional region 
(see Fig. (In this paper we mention several times 
"the inner edge of the peak of the effective potential" ; 
we define this inner edge to be the radius at which the 
effective potential drops to one percent of its maximum 
value at the peak.) 
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FIG. 1: The standing-wave scalar field in a Schwarzschild 
background (solid curve) and the effective potential (dashed 
curve) for angular frequency u> — 0.049. 

Since we are approaching the problem perturbatively, 
we are interested in some small metric perturbation h a p 
on top of the background metric g^ g of Eq. (|14|l that 
would yield the curvature corresponding to the stress- 
energy tensor of the scalar field: 

9a0 = 9% + Kp . (18) 
Linearizing in h a p, this metric gives the Ricci tensor 



R 



a/3 



Ra0 



h V 



1 
2 

and 



x/3\n 



l \af3 



represents the covariant derivative 

,b 

a/3- 



where h 

in the background metric g^ R . For the Schwarzschild 



background metric, R^g = 0. 

We are interested only in spherically symmetric pertur- 
bations. A gauge transformation brings additional sim- 
plification, so h a /3 can be written in the following simple 



/ (l-2/r)H (t,r) 



V o 





H 2 (t,r) 
1-2/r 






0\ 





0/ 



(20) 



(Compare with Eq. (13) of ^1 f° r the case L = 0.) 

We can now substitute h af3 given by Eq. lt2TJ|) into 
Eq. ijT$|) to compute the perturbed Ricci tensor: 



R - ld \ 



2r 3 



d 



d 



{2r~l)—H + —H 2 +r{r 







di 



dr 2 



Rt 



- rdl H2 



d 2 



Ho 



2(r-2) 2 dt 2 2r(r-2) 

3— H Q + (2r-3)— H 2 + r(r 
or or 



Ho 



-2d 
~dr 



Ho 



2 d 



2 dr 



H:> 



(21a) 
#o 
(21b) 
(21c) 
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(21d) 



Inserting expressions (|21|l for R a p into the Einstein 
equations 10, one obtains a set of rather complicated 
PDE's containing both spatial and time derivatives to 
the second order. However, we expect that the equa- 
tions can be further simplified because of additional con- 
sistency conditions imposed on <i> by the wave equation 
(|15fl . Indeed, after adding the R u and R rr equations with 
appropriate coefficients to remove the second derivatives 
in both t and r, and using Rgg = to relate Hq to H 2 , 
we obtain the following set of first-order ODE's for Hq 
and H 2 : 



dH 2 
dr 



Hn 



(r-2y 



dH 
dr 



dH 2 
dr 



-H 



2 ■ 



(22a) 



(22b) 



These far simpler equations can be shown to produce 
no spurious solutions; in fact, together with the wave 
equation (|15|) . they are equivalent to the second-order 
PDE system (|2ll 



2. First-order metric perturbations due to the 
standing-wave scalar field 

In the scalar-field ansatz (|13fl we assumed $ oc coswt. 
Therefore, the driving term on the right hand side of 
Eq. I|22a|) will have static components as well as compo- 
nents oscillating in time at the frequency 2u. Because 
there is no mixing of terms with distinct time signatures 
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m 



equations i|22|) . these terms may be treated separately: 



H 2 (t,r) 



Tjstat 



(r) 



iJ 2 cos (r) cos 2ut; 



H (t,r) = H s tat {r) + H™ s (r) cos 2wf. 



(23a) 
(23b) 



(There is no sin 2wi term with our particular choice of 
the scalar-field phase.) 

For r > 2 analytical approximations for Hq and H% 
are easy to obtain because the scalar field is particularly 
simple there: 



$ w (* /^) cos (tor*) cos (wi), 



(24a) 



where ^0 is the scalar- field amplitude as r — > 00. In- 
serting this into Eqs. I|22|) . we readily compute, at large 




H° tat {r) 

m os {r) 

H^ tat (r) 
H c os (r) 



1 2 T 2 

2 4r 2 



cos 2u>r* 



4r 2 



* 2 1 cos2wr* 
47-2 4r 2 
^wsin 2wr* 

9t9i v]/2 C0S 2 W r* 
^ 2 *olnr + ^^ i 

^gwsin 2u>r* 

47 



'Jjg cos 2wr* 
4r 2 



(24b) 
(24c) 

(24d) 
(24e) 



The static components of Hi and Hq are non-vanishing 
at infinity, and HQ tat actually diverges. This indicates 
that, due to the energy contained in the scalar field, the 
spacetime is not asymptotically flat. However, this bad 
behavior at infinity is an artifact of our model problem 
and is irrelevant to the issues we are studying in this 
paper. 

We can read off from Eqs. (|24|l the ratios of the oscil- 
latory and static components of the metric perturbations 
at large r. They are 



JJCOS 



jjstat 



and 



JJCOS 



Tjstat 



1 

2uor 



Aujt In j 



(25a) 



(25b) 



thus, the static components dominate far from the hori- 
zon. 

Equations l|24|) can be used to set initial conditions for 
the metric perturbations at some large r, allowing for a 
numerical solution to Eqs. ill' 21) from there down to the 
horizon, r — 2. The resulting solution, plotted in Fig. El 
indicates that static components continue to dominate 
near the horizon. 

Near the horizon (inside the effective-potential peak), 
the scalar field has the form 



FIG. 2: Metric perturbations for a standing-wave scalar 
field in a Schwarzschild background with angular frequency 
ui — 0.19 and amplitude = 0.015 far from the black hole, 
corresponding to a binary separation a w 6A/, [Eq. I12al l. 



where "J/^ is the scalar-field amplitude as r — ► 2. Insert- 
ing this approximation into Eq. I|22a|l and averaging the 
right-hand side leads to the following rough estimate of 
the magnitude of the perturbation near the horizon: 



j_j-stat 



2w 2 ^ 2 In (r - 2) 



(26b) 



$ ~ {^in/2) COSOJT* cos cot, 



(26a) 



Inverting this formula can give a useful estimate of the 
distance from the horizon where the perturbation reaches 
a particular value; the estimate turns out to be accurate 
to within a factor of two. 

Although it appears that the metric perturbation di- 
verges at the expected location of the horizon, our per- 
turbative solution is not trustworthy in this regime for 
several reasons in addition to the obvious one of violat- 
ing the perturbative assumption Hq, Hi <C 1: 

1. We ignored the back reaction, i.e., the feedback 
of the metric perturbation into the wave equation. Us- 
ing the Schwarzschild metric in place of the more ac- 
curate perturbed metric in the wave equation, that is, 
using the approximate Eq. (|15|) in place of the exact 
Eq. CP), is equivalent to an error in the scalar- field fre- 
quency Auj/lu « 0(H), which produces phase offsets in 
the scalar field when the wave equation is integrated nu- 
merically. 

2. We linearized the Ricci tensor in the perturbations, 
neglecting higher-order 0(H 2 ) effects. In contrast to the 
linearized equations (|22|) for H2 and Hq, the nonlinear 
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perturbative equations are: 

dH 2 H 2 (l + H 2 ) 

dr r — 2 



(27a) 



(r-2) 2 1-H 



($, t ) 2 + r(l + iJ 2 )($ r ) 2 



dHo 
dr 



= (1 -Ho) 



1 dH 2 



2H 2 



1 + H 2 dr r-2 



(27b) 



Linearization introduces local errors of order H into the 
Einstein equations. However, the errors can build up 
globally when the equations are integrated to obtain a 
numerical solution. The errors produced by linearizing 
the Ricci tensor (the differences between solutions to the 
linearized and nonlinear Einstein equations without back 
reaction in the wave equation) have the same order of 
magnitude in the parameter range of interest as the er- 
rors produced by neglecting back reaction (the differences 
between solutions to the nonlinear Einstein equations de- 
pending on whether wave equation (|f 5fl or Q is used). 

3. We ignored higher harmonics of the scalar field 
and of the metric perturbations that would arise from 
the back reaction. However, these higher harmonics are 
suppressed by additional factors of H cx "J 2 : whereas 
the static and cos 2cut components of H are quadratic in 
^P, higher-order harmonics of frequency 2nu are propor- 
tional to \1/ 2 ™ for n > 1. 



B. Time-averaged fully nonlinear standing- wave 
solution 



To explore the behavior of the standing-wave spheri- 
cal scalar field and the spherical metric near and inside 
the expected location of the horizon, we solve the fully 
nonlinear coupled Einstein-scalar-field equations includ- 
ing full back reaction in the wave equation. To simplify 
our solution, we average the stress-energy tensor in time 
to remove oscillations of the scalar-field energy, so that 
the metric is static. This is justified by the perturbative 
analysis above, which demonstrates that metric compo- 
nents oscillating in time are smaller than static metric 
components and largely decouple from them. 



1. Formalism for nonlinear solution with back reaction 

We write the static spherically symmetric metric in the 
form 

ds 2 = ~e^ +a{r Ut 2 + e f3 ^- a ^dr 2 + rW, (28) 

and we compute the Einstein tensor from this metric in 
the standard way. It is diagonal and its angular compo- 
nents Gqq and are not particularly interesting be- 
cause of spherical symmetry (the angular components of 
the Einstein equation will merely repeat the time and ra- 
dial components by virtue of the contracted Bianchi iden- 
tities). The careted subscripts ^ denote the orthonormal 



basis associated with the (t, r, 9, <f) coordinate system. 
The relevant non-vanishing terms of the Einstein tensor 
in the orthonormal basis are: 



Gff 



z a - p {!3' -a')/r+{l-e a -P)/r 2 , (29a) 
G ff = e a - {P' + a')/r- (1 - e Q ^)/r 2 , (29b) 

where ' denotes a derivative with respect to r, not r* . 

Substituting the Einstein tensor l|29|l and the stress- 
energy tensor (JSJ into the Einstein equations (@J, we ob- 
tain: 



1 



r 







re 



-2a 



K<M 2 



(30a) 



(30b) 



We can now insert the standing-wave scalar-field 
ansatz ff 1 31) and time average the right hand side of 
Eq. (|30b|l over a complete period. For numerical anal- 
ysis it will be useful to switch to a logarithmic coordi- 
nate that changes more gradually than r in the vicinity 
of the Schwarzschild horizon. The following generaliza- 
tion of the Regge- Wheeler tortoise coordinate r* proves 
convenient: 



dr 
dr* 



(31a) 



In terms of this coordinate, the wave equation Q sim- 
plifies to 



d 2 * 2t e a da T 

— — = -uj 2 ^ H — * 

dr* r dr* 



(31b) 



and the Einstein equations 130|) with time-averaged 
(<I>^) 2 and (<& ir ) 2 become 



da 
dr* 



d(3 
dr* 



2r 



* 2 w 2 



dr* 



(31c) 



f 2 e Q * dt> ,„„ ,, 

~ji~nz • (3id) 



2r 3 



r 2 dr* 



2. Singular standing-wave spacetime 

We have solved the coupled equations (|31(l numerically 
to high accuracy for values of the scalar-field amplitude 
and frequency in the range relevant to the BBH problem 
[Eqs. (|12|) ]. Our numerical solutions are very well approx- 
imated by analytic formulae that rely on dividing space 
< r < oo into three regions. Region I is "perturbed 
Schwarzschild", i.e., the region where the perturbative 
solution is valid (r > 2, H < 0.1). Region III describes 
the space very close to r — where the 1/r terms di- 
verge. Finally, the intermediate region II extends from 
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the inner boundary of region I to the outer boundary of 
region III. 

For sufficiently small amplitudes of the scalar field, 
the contributions from the back reaction (by which we 
mean the impact of the deviation of the spacetime from 
Schwarzschild on the solution to the wave equation) and 
from nonlinearity remain small until very close to r = 2, 
so that the metric can be well approximated by perturba- 
tions on top of the Schwarzschild metric. In other words, 
the perturbative solution developed in Sec. 1111 Al is valid 
throughout region I. Indeed, for scalar-field amplitudes 
and frequencies in the range of interest, the metric per- 
turbations HQ tat and if|* at derived in the previous sub- 
section match the values of Hq and Hi corresponding 
to the complete nonlinear solution with back reaction to 
within 3% for H w 0.01. 

We begin the analysis in region III, where r <C 1, by 
assuming e@~ a <C 1 as r — ► 0, which corresponds to 
9rr — * at r = 0. (This assumption, which can be 
deduced from the behavior of d(3/dr* in the transition 
region, will be shown to be self-consistent; more impor- 
tantly, it is supported by our numerical solutions.) Then, 
from Eq. i|31c|) . a' = da/dr — > — 1/r, so a is given by 



a = — In r + ao . 



(32a) 



Here ao is a constant whose value depends on the am- 
plitude and the frequency of the scalar waves; it can be 
roughly approximated by 



a ~ In (A 2 lo 2 ) , 



(32b) 



where A is the amplitude of the scalar field near r = 2. 
The wave equation (|31b(l becomes 

= -^uj 2 e- 2a -a'(W' -tf/r) (32c) 
= -^uj 2 e- 2ao r 2 + l/r(#' - */r) . 

Since we are interested in the region r — > 0, the last term 
dominates, so that the approximate solution to Eq. (|32c() 



"J = nr + kr In r , 



(32d) 



where n, k are constants. 

Substituting \& and a into Eq. (|31d|) and selecting non- 
vanishing terms with the highest order in 1/r, we find 
that ft -> k 2 /(2r), so 



(3= ylnr + /3 , 



(32e) 



where 0o is a constant. Thus, we see that our assumption, 
e P-a ^ ^ as r _> o, is self-consistent: 

(k 2 \ 

(3 — a = f — + 1 j In r + /?q — ao — ► — oo as r — ► 0, (32f ) 

since the coefficient of In r is always positive. 

Our numerical solution in region III agrees well with 
the asymptotic analytical solution (|32|l . For instance, 



the value of k obtained from matching ^ to the form 
of Eq. (|32d() agrees with the value of k obtained from 
matching j3 to Eq. (j32e[) to one part in ten thousand. 
Of particular interest are the metric components and the 
Ricci scalar, whose asymptotics for r — > are: 



9u = 



-a„ r k 2 /2- 



ao r k "72+1 



and 



R 



k 2 e a -l3o r -3-k 2 /2 



(32g) 



(32h) 



(32i) 



The exponent of r in Eq. (|32i() is always negative, so the 
Ricci curvature scalar tends to infinity as r — > 0, i.e., the 
radius of curvature vanishes at the singularity at r = 0, 
as expected. The exponent of r in Eq. (|32h(l is always 
positive, so g rr tends to zero as r — > according to a 
power law. However, the sign of the exponent of r in 
Eq. )32g| ) depends on the value of k, which in turn is 
a complicated function of the scalar-field frequency and 
amplitude. For some scalar field parameter values in the 
range relevant to the BBH problem [Sec.[H] k 2 /2 > 1 and 
gu vanishes at the singularity; for others, gu is infinite 
at r = 0. 

The nature of region II, which represents the transition 
from the Schwarzschild-like region I to the singularity of 
region III, depends strongly on the values of Vfo and co- 
in Schwarzschild, a = ln(l — 2/r) tends to — oo as r — > 2, 
and this is the behavior of a in the nearly Schwarzschild 
region I; meanwhile, in region II, as in region III, a is 
well approximated by 



a 



In r + ao . 



(33a) 



The outer boundary of region II is located at the transi- 
tion between these two behaviors of a, i.e., at the mini- 
mum of a. 

Substituting the approximation l|33a|l for a into the 
wave equation (|31b(l . we obtain: 



d 2 V 
dr* 2 



= y(-LU 2 - 



(33b) 



Thus, the condition for the scalar field to exhibit spatial 
oscillations at an approximately constant amplitude is 
e 2a ° /r 4 <C to 2 ■ The location where this condition begins 
to be violated forms the inner boundary of region II. 
Thus, region II can be said to be defined by the variation 
of a according to Eq. I|33a|l as in region III, and by rapid 
spatial oscillations of the scalar field d^/dr* = to as in 
region I. 

Since ao will be more negative for smaller amplitudes 
of the scalar field, we see that region II is going to be 
significant for small ^o, including those in the relevant 
range of the BBH problem. For larger values of ^0i the 
metric and scalar field will proceed directly from region 
I to region III. 
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When region II does exist, the amplitude and phase of 
the scalar field [solution of Eq. (|33b|l ] 



\Jr(r) = Air) cos </>(r) 



will be given by 



A = A (1 



4r 4 c 



+ -): 



(33d) 



r > 0, albeit with a very large redshift for sources close 
to or inside the location (r = 2) of the Schwarzschild 
horizon. A blown-up view of the region r <C 1 [Figure 



(33c) 4(b) shows that the signal emitted near the singularity 



may be infinitely redshifted or blueshifted depending on 
the asymptotics of the scalar field as r — > according to 



/^rcc e (t»o-/3o)/2 7 .-fc 2 /4+l/2 



f 



(36) 



<i> = u(i + 



2r 4 Lu : 



+ •••), 



(33e) 



to first order in e 2a ° / (r 4 iv 2 ) 

Substituting expressions i|33|) for a and \& into the dif- 
ferential equation for (3, Eq. (|31d(l . we find that the dom- 
inant term is the first one, dfl/dr* 
in region II (3 is approximately 

f 



(l/2)e- a °A 2 w 2 



so 



(3 = ^e- a °A 2 Lu 2 r* + const 



~ 2ao r 2 A 2 uj 2 + const. 



(33f) 



where the last equality comes from the integral of equa- 
tion ifSTajl . r* = e~ Q °r 2 /2 + const. 

Embedding diagrams and redshifts may provide 
the best pictorial insight into our full time-averaged 
standing-wave scalar-field solution, including all of re- 
gions I, II and III. 

Figure shows an embedding diagram for the 
standing-wave spacetime: 



-T- = V\9rr - 1| 



(34) 



The 2-surface obtained by rotating around the vertical 
axis r = has the same 2-geometry as the surface (t — 
const, 6 = 7r/2) in the standing- wave spacetime. At radii 
r > 2 the embedding is very nearly the same as for a 
Schwarzschild black hole (cf. Fig. 31.5 of 7]). For r < 2, 
the radial distance changes far more slowly than the areal 
radius (0 < g rr <C 1), so the embedding is performed in 
Minkowski space rather than Euclidean space: the metric 
is ds 2 — —dz 2 + dr 2 + r 2 d(j) 2 rather than ds 2 — +dz 2 + 
dr 2 + r 2 d(f> 2 . The embedded surface asymptotes to the 
light cone as r — > 0. 

Figure 0] depicts the redshift of light emitted at one 
radius and received at another, greater one, as a function 
of the emitting radius: 



(35) 




Figure 4(a) shows that, while the redshift becomes very 
large as r — > 2, it never becomes infinite there as it would 
for a Schwarzschild black hole. As expected, the horizon 
is destroyed by the standing-wave scalar field, so an ob- 
server at infinity can receive signals from any source at 



in .-■ 



9 10 



FIG. 3: Embedding diagram for the spacetime with time- 
averaged standing-wave scalar field of angular frequency u) = 
0.19 and amplitude "Jo = 0.015 at large radii [correspond- 
ing to the binary black hole separation a w 6M; Eq. 112aH l. 
The solid line represents embedding in Euclidean space; the 
dashed line, embedding in Minkowski space. Regions I, II and 
III are labeled on plot. 



3. Comparison of standing-wave and Schwarzschild 
spacetimes 

It is important to understand how the complete 
standing-wave spacetime with back reaction (we shall call 
this spacetime S) compares with the Schwarzschild space- 
time (which we shall call spacetime D). We might first 
try to compare the metric components in the two space- 
times as functions of the radial coordinate r. Indeed, 

" are 



r 2 sin 2 i 



the metric components gee — r 2 and 
precisely equal in the two spacetimes when evaluated at 
the same location in (t, r, 9, <f>) coordinates. Furthermore, 
outside the effective-potential region, the perturbation 
due to the scalar field is so small that the fractional dif- 
ference 5g a p/g al3 = (g%p - g^)/g^ in metric compo- 
nents gtt and g rr does not exceed 0.01% for scalar-field 
parameters in the range of interest. However, the metric 
components gtt and g rr in S and D can differ by orders 
of magnitude near r = 2, inside the effective potential 
peak. 

The apparent mismatch between the metric compo- 
nents of the two spacetimes near r — 2 turns out to be 
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FIG. 4: (a) Redshift z = 8X/X of light emitted from radius 
r and received by an observer at r = 10. (b) Redshift for 
an observer at r = 0.0001. A distant observer would see 
light emitted from r = 0.0001 redshifted by \n(z + 1) « 10 5 . 
These curves are drawn for the spacetime with time-averaged 
standing-wave scalar field that has angular frequency 10 = 
0.19 and amplitude *l/o = 0.015 at large radii [corresponding 
to the binary black hole separation a ~ 6M; Eq. H12al 1. 



a consequence of a poor choice of the radial coordinate 
r for comparison. A much better choice is r*\ when the 
coordinates (t, r* , 0, <j>) are used for mapping between the 
two spacetimes S and D, the metric components agree 
extremely well near r = 2. 

The fractional differences Sg/g between the gtt and gee 
components in S and D are plotted in Fig. [5] for scalar 
field parameters corresponding to binary black hole sep- 
arations at the boundaries of the range of interest. Using 
r* rather than r as the coordinate for comparison means 
that the gee components no longer match perfectly; how- 
ever, the fractional difference introduced remains small as 
r — > 2 and does not exceed 0.6% in the range of interest. 
The fractional differences in g^ are identical to those in 
gee and are not plotted separately. The Regge- Wheeler 



tortoise coordinate r* [Eq. 1|16[) ] and its generalization 
[Eq. Q31a[l ] were defined so that g r * r * = —gtt in both 
spacetimes S and D; therefore, the fractional differences 
in the values of g r * r * in S and D are the same as the 
fractional differences in g tt . 

As Fig. [5] shows, the fractional differences in the met- 
rics are < 0.02 down to values of r* ~ —1000, a location 
so deep inside the peak of the effective potential that it 
contains at least 20 near-horizon oscillations of the scalar 
field for frequencies and amplitudes in the BBH separa- 
tion range of interest. Perhaps a more impressive way 
to state this is that in the (t, r* , 8, <j>) coordinate system, 
metric components of g s and g D match to an accuracy 
of 2% for all relevant scalar-field parameters down to the 
Schwarzschild radius r D — 2 < 10 -100 . 

The fractional differences between the coefficients of 
the metrics g s and g D continue to grow approximately 
linearly in r* deep inside the effective potential and reach 
10% at the Schwarzschild radius r D ~ 2 - 10~ 3000 , 
or approximately 500 scalar-field oscillations inside the 
effective-potential peak for scalar field parameters corre- 
sponding to BBH separation a w 6M. 



IV. DOWNGOING SCALAR FIELD 



Having discussed, in Sec. lIIII the standing- wave scalar- 
field spacetime that modeled the stationary BBH approx- 
imation, we now turn to a scalar-field spacetime that 
serves as a model for the physical BBH spacetime with 
downgoing gravitational waves at the black-hole hori- 
zons: a nearly Schwarzschild spacetime with spherically 
symmetric scalar waves that are purely downgoing at 
r = 2. 

For a perturbative analysis of downgoing scalar waves 
in Schwarzschild, the ingoing, null Eddington-Finkelstein 
time coordinate v = t + r* is more appropriate than 
the standard Schwarzschild time coordinate. Let us 
suppose that by the time v = a total mass-energy 
Mq = m(v = 0) is located within the horizon r = 2. 
We are not particularly interested in how this mass ac- 
cumulated there or how the scalar field behaved in the 
past; we are only interested in the times immediately 
following v = 0, and we let the scalar waves be purely 
downgoing and monochromatic at the horizon for v > 0. 
Then for v > radiation is falling into the black hole 
at a nearly constant rate, corresponding to the energy 
density in the scalar field dm/dv w \E , go; 2 /2, with some 
small oscillations on top of the linear increase in mass. 
This is very similar to the Vaidya solution and, indeed, 
the Vaidya metric will be seen to describe the spacetime 
of the downgoing scalar-field solution: 



ds 2 = - 



1 - 



2m(v) 



dv 2 + 2dvdr + r 2 dtt 2 . 



(37) 



Near the horizon, $ = (1/r) cos luv is a purely down- 
going solution to the wave equation |{T{|, The only non- 
zero term of the Ricci tensor in Vaidya coordinates is 



10 



0.02 



0.015 



0.01 - 



0.005 - 



1,1,1,1,1,1 

















500 -250 



250 
r* 



500 



750 



1000 



the Schwarzschild effective-potential peak. Outside the 
effective-potential peak there is both a downgoing scalar 
field and an upgoing one, reflected off the potential. Since 
for small v the metric is nearly Schwarzschild [the con- 
stant Schwarzschild mass M is replaced by the m(v) 
of Eq. (32} ], the scalar field everywhere is given to a 
high accuracy by a solution to the wave equation in the 
Schwarzschild background subject to the purely downgo- 
ing boundary condition at the horizon. (Of course, very 
far from the horizon the energy contained in the inter- 
vening scalar field will act as an additional mass, but we 
are not interested in this region for our model problem.) 



V. RECONSTRUCTION OF DOWNGOING 
SCALAR FIELD FROM STANDING- WAVE 
SCALAR FIELD 
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FIG. 5: (a)Fractional differences of the metric components 
gu ~ —g r 'r' (solid curve) and gee (dashed curve) between 
Schwarzschild spacetime D and standing-wave scalar field 
spacetime S with scalar-wave amplitude and frequency cho- 
sen to model BBH separation a « 6M [Eq. 112al l. (b)Same 
quantities plotted for scalar field parameters chosen to model 
BBH separation a w 15M [Eq. iT2bl l. 



R vv — (2/r 2 )m'(v), where ' denotes the derivative with 
respect to v. The Einstein equations © at r = 2 say: 



2m' (u) 



2^qUJ 2 sin ojv 



4 •» - 4 

Equation (|38Jl is trivially integrated to obtain: 

^qUJ sin 2ojv 



m(v) = M + 



(38) 



(39) 



The black-hole mass grows linearly in v at the rate 
^uj 2 /2 with a tiny superimposed oscillatory component. 
The black hole retains a smooth, non-singular future 
horizon. 

The scalar field is purely downgoing at the hori- 
zon and approximately downgoing everywhere inside 



We turn now to our scalar-wave version of adding a 
radiation-reaction field to a standing-wave spacetime to 
obtain a physical spacetime with downgoing waves at 
horizons and outgoing waves at infinity. For this pro- 
cedure there is a substantial difference between the BBH 
problem and our model problem. 

In the true BBH problem, the periodic standing wave 
(SW) solution is sourced by the black holes and corre- 
sponds to the ^Retarded + ^Advanced solution of the 
Green's function problem. In this case we add the non- 
sourced ^Retarded — ^Advanced radiation reaction (RR) 
solution of the linearized Einstein equations in the SW 
spacetime to get an approximation to the physical re- 
tarded solution At infinity, where the SW field is 
^Outgoing + ^Ingoing, the boundary condition for the 
RR field should be set to \ Outgoing — ~ Ingoing, so 
that their sum contains only physical outgoing waves, 
and similarly at the horizons the RR field will be 
^Downgoing — ^Upgoing. Adding this RR field to the 
^Downgoing + ^Upgoing standing waves would yield 
gravitational waves that are downgoing at the expected 
horizon locations, conforming to the expected behavior 
in physical black-hole spacetimes. (We do not expect 
the stress-energy tensor of the sum of SW and RR waves 
to precisely match the Einstein tensor of the SW space- 
time because, of course, gravitational theory is not linear; 
however, it is likely that "effective linearity" holds in the 
sense defined by Price for the non-tidally-locked case 
as well as for the tidally-locked case. In a future paper 
we intend to explore this issue with a model that more 
closely resembles the BBH problem.) 

The scalar-field model we are currently analyzing is 
not sourced: the wave equation Q we used to compute 
the SW solution is homogeneous. There is then no per- 
turbative homogeneous solution that is simultaneously 
^Outgoing — ^Ingoing at infinity and ^Downgoing — 
^Upgoing at the expected horizon location. Since at the 
outer boundary the problem is obviously linear for suf- 
ficiently weak scalar fields, it is easy to reconstruct the 
outgoing solution from the SW solution there: we sim- 
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ply double the outgoing component of the SW solution. 
The interesting case lies in the extraction of a downgo- 
ing solution near r ~ 2. We attempt to reconstruct the 
downgoing scalar field from the SW scalar field near the 
expected horizon by adding to the SW field a perturba- 
tive RR field that is iDowngoing — i Upgoing at r ~ 2. 



As in Sec. IIII Bl let S denote the spacetime of the com- 
plete standing-wave solution with back reaction. As dis- 
cussed in the previous section, the spacetime of the down- 
going scalar field is approximated to sufficient accuracy 
for our purposes by the Schwarzschild spacetime D. 

The complete SW scalar field is a solution to the wave 
equation in spacetime S (in our simplified treatment of 
the problem, spacetime S actually corresponds to the 
time-averaged solution, i.e., one in which we ignore the 
oscillatory components of the metric). The RR field is a 
solution to the same wave equation in S in our model. 
The "reconstructed" downgoing field is, therefore, the 
downgoing solution to the wave equation in S. We want 
to compare this to the "true" downgoing field, which is 
the downgoing solution to the wave equation in D, i.e., 
in Schwarzschild. 

In the region between the expected horizon location 
r = 2 and the inner edge of the peak of the effective 
potential, the wave equation is dominated by 



d 2 ^ 
dr* 2 



(40) 
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in both spacetimes S and D. Hence, the solution to the 
wave equation will be oscillatory in r* with frequency 
uj, which makes sense on physical grounds, since ingoing 
light cones are t + r* = constant in both S and D. More- 
as discussed in Sec. 



over, as discussed in bee. IIII Bl the metrics of the two 
spacetimes are nearly the same in the r* coordinate, i.e., 
g s (r*) w g D (r*). This suggests that to get the scalar 
wave phasing to agree, we need to map between the two 
spacetimes using the r* radial coordinate. 

We set the boundary conditions for both the RR field 
in S and the downgoing field in D at a negative value of 
r* chosen so that the fields are at least a few wavelengths 
inside the effective potential, and so that r s (r*) is very 
close to r s = 2 (it might actually be slightly inside r = 2). 
The SW + RR and downgoing scalar fields will match by 
construction at the point where the initial conditions are 
set. We will integrate both solutions toward larger r* 
and compare the quality of the match between the two 
fields. 

For the purposes of comparing the scalar fields in the 
two spacetimes, we separate the complex scalar field 
fy(r*) [the spatial factor of the complete field $(r, t) — 
ft[\f r (r*)e~ ia '*]/r, cf. Eq. (J3J] into upgoing and downgoing 
components. We define the amplitudes and phases of the 
upgoing and downgoing fields as follows (see below for 



FIG. 6: (a) The fractional difference in the amplitudes of the 
reconstructed scalar field and downgoing scalar field 5d/d 



{d sw+i 



d" 



l )/d d 



(solid curve) and the phase differ- 



ence between the two fields 5<f> d = <t> S d W+RR - (j>f wn (dashed 
curve), plotted vs. r*. Scalar- wave amplitude and frequency 
chosen to model BBH separation a m 6M. (b) Same quanti- 
ties plotted for scalar-field parameters chosen to model BBH 
separation a m 15M. 



motivation) : 



d = — 



I 
1 
f 



dr* 
dr* 



2icuu V dr* 
1 /d* 



—2iujd \ dr 



(41a) 
(41b) 
(41c) 
(41d) 



To motivate these definitions we consider the geomet- 
ric optics limit, where the wave phase evolves much faster 
than the amplitude. In this limit, the downgoing compo- 



12 



nent of the scalar field ^>d oc e separates unambigu- 
ously from the upgoing component ^> u cx e lwr : 



*(r*) = we 10 " +de^ d , 
where we use the standard approximations 



d(j) u 








= W > 


dr* 


dr* 


d(f>d 








= w > 


dr* 


dr* 



(42a) 

(42b) 
(42c) 



Inverting Eq. I|42a|) with these approximations yields the 
definitions l|41(l . Although the geometric optics approx- 
imations break down in the region of the effective po- 
tential, and the separation of the scalar waves into up- 
going and downgoing components becomes ambiguous 
there because the wave speed is ill-determined outside 
the short- wavelength limit, expressions (|41|1 are adequate 
for comparing scalar fields in our region of interest. 

In Fig. we show the fractional difference Sd/d = 
( d SW+RR _ d downy d down in the amp Htude of the down- 
going components of the reconstructed SW + RR waves 
and the downgoing waves along with the phase differ- 
ence 5(f>d = 4>l W+RR - (f>j own - The two plots repre- 
sent the endpoints of t he ran ge of relevant BBH sepa- 
rations: a 6M in Fig. 6(a) and a w 15M in Fig. 6(b) 



Only the downgoing amplitude d and downgoing phase 
4>d are plotted. The upgoing field components are zero 
to numerical precision inside the effective potential and 
the differences between the reflected upgoing components 
of the "reconstructed" and "true" downgoing fields out- 
side the effective-potential peak are similar to the dif- 
ferences between the downgoing field components there, 
Su/u ~ Sd/d and 8(j) u ~ S(pd- 



The amplitudes and phases of the "true" downgoing 
field and the "reconstructed" downgoing field match to 
within one part in ten million from the location where 
the initial conditions are set (several scalar-field oscilla- 
tions inside the effective potential) to the inner edge of 
the effective-potential peak for all BBH separations in 
the range of interest. Near the effective-potential peak 
the fractional difference in the amplitudes does not ex- 
ceed 0.03% and the phase difference is less than 0.002. 
Outside the effective potential, the fractional difference 
in the amplitudes is 5 parts per million and the phase 
difference is less than 0.00002 for the smallest BBH sep- 
arations in the range of interest. 

We also compared the "reconstructed" and "true" 
downgoing fields very deep inside the effective potential 
when the field-matching initial conditions are set about 
10 scalar-field oscillations inside the effective-potential 
peak. In this case, the amplitudes of the two fields are 
equal to within numerical precision and the phase differ- 
ence does not exceed 3 x 10~ 7 down to 500 scalar-field 
oscillations inside the effective-potential peak. The fields 
begin to disagree significantly only once the naked singu- 
larity is approached in the spacetime S, at r s (r*) < 0.2 



0. 
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